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1. INTRODUCTION 

Interest in graphene-based superlattices has increa- 
sed in recent years. Calculations of graphene-based 
superlattices with periodic rows of vacancies were per- 
formed using the molecular dynamics method [1] . Cal- 
culations of single-atom-thick superlattices formed by 
lines of pairs of adsorbed hydrogen atoms on graphene 
were carried out with the density functional theory [2] . 

Rippled graphene that can be treated as a superlat- 
tice with the one-dimensional periodic potential of rip- 
ples was investigated in [3-5]. Superlattices obtained 
when a periodic electrostatic potential [6-9] or periodi- 
cally located magnetic barriers [10-13] were applied to 
graphene were analytically examined. 

However, the investigation of the graphene-based 
superlattice with a periodic electrostatic potential dis- 
regarded the fact that the application of the electro- 
static potential to a gapless semiconductor (graphene) 
results in the production of electron-hole pairs and the 
redistribution of charge s: electrons move from the 
region where the top of the valence band lies above 
the Fermi level to the region where the bottom of the 
conduction band lies below the Fermi level. The su- 
perlattice becomes a structure consisting of positively 
charged regions, where the electrostatic potential dis- 
placing the Dirac points upward in energy is applied, 
alternating with negatively charged regions. The strong 
electrostatic potential of induced charges appears and 
strongly distorts the initial step electrostatic potential 
and, therefore, the electronic structure of the superlat- 
tice calculated disregarding the electrostatic potential 
of induced charges. 

To avoid the production of electron-hole pairs, a 
superlattice appearing due to the periodic modulation 
of the band gap is considered. 



A superlattice in the form of the periodic planar 
heterostructure of graphene nanoribbons between which 
nanoribbons of hexagonal boron nitride (h-BN) are in- 
serted was previously proposed in [14]. The band struc- 
ture of such a superlattice was numerically calculated. 
However, it is very difficult to implement this super- 
lattice even using the advances of modern lithography, 
because problems inevitably arise with the control of 
periodicity in the process of the etching of nanoribbons 
in a graphene sheet and the insertion of h-BN nanorib- 
bons. Moreover, h-BN is an insulator with a band gap 
of 5.97 eV, which significantly hinders the tunneling of 
carriers between graphene nanoribbons. Such a het- 
erostructure is most probably a set of quantum wells 
where the wavefunctions of carriers from neighboring 
quantum wells almost do not overlap. 

In this work, a superlattice formed by a graphene 
sheet deposited on a strip substrate is proposed. The 
strip substrate is made of periodically alternating strips 
of Si02 (or any other material that does not affect the 
band structure of graphene) and h-BN, as shown in 
Fig. 1. The h-BN layers are located so that its hexag- 
onal crystal lattice is under the hexagonal crystal lat- 
tice of graphene. Owing to this location, a band gap 
of 53 meV appears in the band structure of graphene 




Fig. 1. Graphene on the strip substrate consisting of 
alternating Si02 and h-BN strips. 



Fig. 2. One-dimensional periodic Kronig-Penney 
potential of the supciiatticc show n in Fig. 1: the 
periodically alternating gap modification of graphene 
on h-BN with a band gap of 2Ao = 53 meV and 
gapless graphene on SiC>2. 

in the graphene-sheet regions under the h-BN layers 
[15, 16] (graphene with a band gap is called a gap mod- 
ification of graphene) . 

It is assumed that all of the contacts between the 
regions of different band gaps are first-kind contacts 
(Dirac points of graphene are located in the band gaps 
of the gap modification of graphene). Such a superlat- 
tice is a first-type superlattice (classification of super- 
lattices can be found, e.g., in [17]). 

There are other gap modifications of graphene. A 
band gap in graphene epitaxially grown on the SiC sub- 
strate is nonzero [18] and, according to experiments 
with the use of angular resolved photoemission spec- 
troscopy, is 0.26 eV [19]. An additional modification of 
graphene, graphane [20], where the direct band gap at 
the point T is 5.4 eV according to the calculations [21], 
was recently synthesized by hydration. However, the 
insertion of epitaxially grown graphene between gap- 
less graphene strips is not easier than the insertion of 
h-BN. The band gap of graphane is too wide. These 
gap modifications of graphene are inapplicable for su- 
perlattices. 

The main advantage of the proposed superlattice 
is the simplicity of the manufacture and control of 
its periodicity. It is worth noting that some prob- 
lems can arise in the superlattice. The difference be- 
tween the lattice constants of h-BN and graphene is 
about 2% [16]. If about 100 hexagonal graphene cells 
are packed into one period of the superlattice, the for- 
mation of the band gap in the gap modification of 
graphene in the graphene sheet regions above h-BN is 
violated owing to the inaccurate arrangement of carbon 
atoms above boron or nitrogen atoms. Since contacts 
between graphene and its gap modification are not het- 
erocontacts (contacts between substances with differ- 
ent chemical compositions) , the edges of quantum wells 
can be insufficiently sharp and quantum wells cannot 
be considered as square quantum wells. A transient 
layer with a spatially varying band gap can exist in- 
stead of the sharp edge. Finally, the substrate can be 
stressed. The appearing periodic stress field of the sub- 



strate can also affect the band structure of the proposed 
superlattice, but this effect is very small. 

2. MODEL DESCRIBING THE SUPERLATTICE 

The x and y axes are perpendicular and parallel 
to the interfaces of h-BN and Si02 strips, respectively 
(see Fig. 1). The superlattice is described by the Dirac 
equation 

{v F( rp + Aa z + V)t>(xy) = E^{x, y), (1) 

where vf ~ 10 8 cm/s is the Fermi velocity, cr = (a x , a y ) 
and <j z are the Pauli matrices, and p = — iV is the 
momentum operator (the system of units with h = 1 is 
used) . The half- width of the band gap is periodically 
modulated: 



A = 



0, d(n — 1) < x < —djj + dn, 
Aq, — dn + dn < x < dn, 



where n is an integer enumerating the supercells of the 
superlattice; di and dn are the widths of the Si02 and 
h-BN strips, respectively; and d = dj + d H are the pe- 
riod of the superlattice (the size of the supercell along 
the x axis). The periodic scalar potential V can appear 
due to the difference between the energy positions of 
the middle of the band gap of the gap modification 
of graphene and conic points of the Brillouin zone of 
gapless graphene (see Fig. 2): 



V = 



0, d(n — 1) < x < — djj + dn, 
Vq, —da + dn < x < dn. 



In order for the superlattice to be a first-type super- 
lattice, the inequality |Vo| < should be satisfied. The 
solution of Eq. (1) for the first supercell has the form 

*(x, y) = ipx(x)e ikvV , < x < d. 

For the nth supercell, in view of the periodicity of the 
superlattice, 

ipn{%) =i*i(x+(n— l)d). 

In the quantum well region (0 < x < dj), the solution 
of Eq. (1) is a plane wave 



V£ 1} W = N kl (jw)e*» + N kl (^)e-*«, (2) 

where N kl is the normalization factor. The substitu- 
tion of Eq. (2) into Eq. (1) provides the relation be- 
tween the lower and upper spinor components 

= A^a« rfW = -A c (1) A+ = v F^ ±ik v) 

The relation of E with k\ and k y has the form 



E = ±v F Jkl + kl 



2 



It is convenient to represent Eq. (2) in a more compact 
form [7]: 



These conditions provide the equalities 



i {2} \ / 



^1 r (i) 



(3) 



(i) 

'n+l 

.(1) 
n+l 



n^(o)n fc2 (d)( 

^ Cn 



When the inequality 



According to definition (7) of the T matrix and the last 
two equalities 1 , 



A 2 + v 2 F k 2 -(E-V ) 2 >0 



(4) 



is satisfied, the solution of Eq. (1) in the barrier region 
(dj < x < d) has the form 



t = n^ 1 (o)n fca (d)n^ 1 (d / )n fcl (d / ). 



(8) 



The substitution of Eqs. (3) and (5) with the corre- 
sponding arguments into Eq. (8) yields the expressions 



^(x) = n k2 (x) r? 2) 



(2) 



n k2 { x ) = N k2 



where 



-A_ A. 



1 



(5) 



ae 



ik\dj 



(A_+A + )(A+ + A_)e- fe2d "- 
-(A_-A_)(A + -A + ) e fe ^»], 
2ae~ lkldl (\- + A+)(A_ - X_)sh(k 2 dn), 



k%x a z 



(9) 



T, 



~ x - = IfS^' fc2 = ^ + « - - ^ — 
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The solution of Eq. (1) in the barrier region under the 
condition 



A 2 +v 2 F k 2 y -(E-V ) 2 <0 



(6) 



is given by Eq. (5) with the change ki —¥ i*t%, i.e., it 
is oscillating. 

The possibility of existing Tamm minibands formed 
by localized states near the interface between graphene 
and its gap modification will be considered below. In 
this case, k\ — > i%\ k% is real. A necessary condition 
for existing Tamm states has the form 

\ k y\ > W Ii 



under this condition, the energy E = ±v F ^J k^ — xf is 
real. 

3. DERIVATION OF THE DISPERSION 
RELATION 

The dispersion relation is derived using the transfer 
matrix (T matrix) method. The T matrix relates the 
spinor components for the nth supercell to the spinor 
components of the solution of the same type for the (n 
+ l)th supercell. For example, for the solution in the 
quantum well region, 



(A + + A_)(A+ + A_) 

The last two relations in Eqs. (9) are the general 
properties of the T matrix. 

The derivation of the dispersion relation with the 
use of the T matrix is briefly as follows. 

Let N — L/d be the number of supercells in the 
entire superlattice, where L is the length of the super- 
lattice along the x axis, i.e., the direction of the ap- 
plication of the periodic potential. The Born-Karman 
cyclic boundary conditions for the superlattice have the 
form 

^ 2 \x)=^' 2 \x). 
At the same time, 

^ 2 \x)=T^ 2 \x), 

from which, T N = I, where I is the 2x2 unit matrix. 

It is convenient to diagonalize the T matrix by 
means of the transition matrix S: 



T rf = STS- 1 = 



Ai 






A 2 



a (1) \ 

"n+l \ _ ( a n 

Ji) J ~M c w 



(7) 



To determine the T matrix, the following conditions 
of the continuity of the solution of the Dirac equation 
describing the considered superlattice are used: 

4P(d I -o) = 4 2 Hd I + o), 
^ 2 Hd-o) = ^l(+o). 



where Ai^ are the eigenvalues of the T matrix and have 
the property A 2 = . According to = I 

Ai = e 2m/N , -N/2 <n< N/2. 

In view of the property TrT = TrTd and in terms 
of the notation k x — 2irn/L (— 7r/d < k x < ir/d), the 
dispersion relation is obtained in the form 



TrT = 2cos(k x d). 



(10) 



1 Note that the cyclic permutations of the factors of Q matrices 
are possible in the definition of the T matrix; these permutation 
do not change dispersion relation (10). This can be verified by 
comparing Eq. (8) with Eq. (23) in [7]. 
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) -sh(k,2dn) sin(feidj) 



(11) 



Taking into account the last relation in Eqs. (9), Eq. (10) 
can also be written in the form 

ReTu — cos(k x d). 

Dispersion relation (10) under condition (4) gives the 
equation 

v\k\ - v\k\ + V 2 - A 2 
2vpkik2 

+ ch(k2du) cos(kidi) = cos(k x d). 

According to this equation, the passage to the single- 
band limit is performed by two methods: hrst, Vq = Aq 
(the quantum well only for electrons) and, second, Vo 
= — A (the quantum well only for holes). The result 
of the passage coincides with the known nonrelativistic 
dispersion relation (see, e.g., [22]), although the expres- 
sions for fci, &2, and E arc different. 

If inequality (6) is satisfied, the change k2 — > 1x2 
should be made in Eq. (11): 



-v 2 F x 2 - vpk\ + V 2 - A 2 a 



sin(>f2rfjj) sin(fcidj) 



2v 2 F kix 2 v ' v ' " (12) 

+ cos(x2du) cos(k\di) = cos(k x d). 

For Tamm minibands, the change ki — > ia\ should be 
made in Eq. (11): 



vpk\ + v 2 F x\ + Vq - A 



2vpXi k2 



- sh{k2di 1) sh(x\di) 



(13) 



+ ch{k2di i)ch{x\di) = cos(k x d). 

Equation (13) has the solution only under the condition 

v 2 F k 2 2 + v 2 F x 2 + V£ - A 2 < 0. 

However, this inequality is never satisfied; correspond- 
ingly, the Tamm minibands are absent in the superlat- 
tice under consideration. 

4. RESULTS OF THE NUMERICAL 
CALCULATION 

The numerical calculations of the dependence of the 
energy on k x were performed for two values k y = and 
0.1 imi 4 at Vo = (see Fig. 3). The energy of car- 
riers is assumed to be low, \E\ < 1 eV, because the 
Dirac dispersion relation for carriers and, correspond- 
ingly, Dirac equation (1) are invalid for high energies. 

The electron minibands are separated from hole 
minibands by a band gap, which increases with \k y \. 
For d[ = djj at k y = 0, it is E g ~ 10-30 meV when 
d=10 — 100 nm. In this case, the solution of Eq. (11) is 
transformed to the solution of Eq. (12). The band gap 
can increase strongly when djj increases with respect 
to dj: Eg > 100 meV, i.e., is several times larger than 
2A . 

The width of the minibands decreases with an in- 
crease in the period of the superlattice d. The depen- 
dence of the width of the minibands on Vq was also 
examined. The widths of the electron and hole mini- 
bands increase and decrease, respectively, at Vo > 
and vice versa at Vq < 0. 
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Fig. 3. Numerically calculated dependence of the 
energy on k x for two k y values and two superlattice 
periods d. The dispersion curves for the superlattices 
with di = (solid lines) du, (dashed lines) du/2, and 
(dotted lines) 2djj. 
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5. POSSIBLE APPLICATIONS 
OF THE SUPERLATTICE 

The described superlattice can be used as a field- 
effect transistor (FET) where the substrate serves as a 
gate. The ratio of the current through the superlattice 
to the current through the gate at a substrate thickness 
of about 10 nm can reach ~ 10 6 as for FET based on 
graphene nanoribbons [23] . The main advantage of the 
considered superlattice is the absence of the effect of 
the scattering of carriers on the edges of a nanoribbon 
on their mobility. The mobility of the carriers in gap- 
less graphene reaches /x = 2 x 10 5 cm 2 /(V s) [24, 25]. 
However, the mobility of carriers in FET based on the 
graphene nanoribbon with a width of w ~ 3 nm is three 
orders of magnitude smaller than fiQ. The cause of 
such a strong decrease is possibly the scattering of car- 
riers at the edges of graphene nanoribbons. The mean 
free path between two acts of the scattering of carriers 
at the edge of the graphene nanoribbon X e d ge w/P, 
where P is the probability of backscattering [23] . For 
sufficiently good edges, P<Cl. The problem of scatter- 
ing on edges is absent for the proposed superlattice; 
for this reason, the mobility of the carriers in the su- 
perlattice is expected to be ~ /i m the absence of the 
problems with the periodicity of the potential. At the 
same time, a sufficiently large E g value, which pro- 
vides the operation of FET at room temperature, can 
be reached. 

If an Au film is deposited on the lower side of the 
substrate and graphene is optically pumped, the su- 
perlattice can be used as a terahertz laser similar to 
a terahertz laser based on gapless graphene [26]. In 
this case, terahertz radiation will be emitted from the 
regions of the Si02 substrate. 



6. CONCLUSIONS 

A model describing such a superlattice based on 
graphene on a strip substrate has been proposed. The 
dispersion relation has been derived, which is trans- 
ferred to the known nonrelativistic dispersion relation 
in the passage to the single-band limit. The numeri- 
cal calculations have been performed for a pair of the 
nearest electron and hole minibands using the derived 
dispersion relation. Possible applications of the super- 
lattice as a transistor or a terahertz laser have been 
pointed out. 
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